Interferometry radii in heavy-ion collisions at y/s = 200GeV and 2.76TeV 
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The expansion of the firebaU created in Au-Au collisions at = 200GeV and Pb-Pb colhsions 
at 2.76TeV is modelled using the relativistic viscous hydrodynamics. The experimentally observed 
interferometry radii are well reproduced. Additional pre-equilibrium flow improves slightly the 
results for the lower energies studied. 



PACS numbers: 25.75.Ld, 24.10Nz, 24.10Pa 

Keywords: relativistic heavy-ion collisions, hydrodynamic model, femtoscopy 



o 



> 

in 



I. INTRODUCTION 

The matter created in ultrarelativistic heavy-ion col- 
lisions undergoes a rapid expansion. During the expan- 
sion substantial transverse flow is generated modifying 
the transverse momentum spectra of emitted particles. 
The measured spectra and elliptic flow at the BNL Rel- 
ativistic Heavy Ion Collider [1] can be described within 
relativistic hydrodynamic models [2, 3]. The hydrody- 
namic behavior of the bulk of the matter created in Pb- 
Pb collisions at ^/s = 2.76TeV is confirmed in recent 
experiments at the CERN Large Hadron Collider (LHC) 
[4, 5]. The magnitude of the observed elliptic flow indi- 
cates that the dense matter behaves as an almost perfect, 
thermalized fluid, with small shear viscosity [6]. 

The size and the life-time of the source can be esti- 
mated from the Hanbury Brown-Twiss (HBT) radii ex- 
tracted from Bose-Einstein correlations of identical par- 
ticles [7] . The description of the HBT radii in dynamical 
models has been the subject of intensive studies. Many 
calculations lead to an overestimate of the ratio of Rout to 
Rside radii, which was termed as the HBT puzzle. Cap- 
turing correctly the values of the Rside and Riong radii in 
a model calculation amounts to a satisfactory description 
of the the size and the life-time of the system. The radius 
Rout (and the ratio Rout / Rside) are more sensitive to the 
collective flow of the matter at the freeze-out. The same 
is true for the dependence of the radii on the momen- 
tum of the pion pair. Experimental observations yield 
Rout /Rside — Ij contrary to expectations for a system 
undergoing a first order phase transition [8]. Reversing 
this argument, the experimental data indicate that the 
equation of state of the hot matter is hard, without a 
soft point [9, 10]. Other effects that influence the ratio 
Rout /Rside are the shear viscosity [10, 11], the partial 
chemical equilibrium [12], the early start of the hydrody- 
namic expansion (0.1-0.25fm/c) [9, 10] or the presence of 
some pre-equilibrium flow [13]. 

We simulate central heavy-ion collisions at the top 
RHIC and LHC energies using relativistic viscous hydro- 



dynamics, with or without pre-equilibrium flow. We cal- 
culate the transverse momentum spectra and HBT radii. 
We obtain a satisfactory agreement with the experiments, 
the additional initial flow improves the Rout /Rside ratio 
at RHIC energies but has a smaller effect for Pb-Pb col- 
lisions at the LHC. 



II. HYDRODYNAMIC MODEL 

We use the second order relativistic viscous hydrody- 
namics [14] to model the expansion of fireball. In rel- 
ativistic viscous hydrodynamics the energy-momentum 
tensor is 



T^"^ = (e + - pg'"' + Ti^"" + HA^ 
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with stress corrections from shear tt'*" and bulk HA^'' 



viscosities, A^" — g^" 



The stress corrections 



are the solutions of the second order relativistic viscous 
hydrodynamic equations [6, 11, 14, 15] 
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Tj and C are the shear and bulk viscosity coefficients 
and two relaxation times r^r and rn appear. We use 
the Navier-Stokes initial conditions for the stress tensor 
27r^^(ro) = 27r^^(To) — — 7r^^(To) = 47^/3x0 corresponding 
to the longitudinal Bjorken flow and H(to) = 0. 

The initial profile of the entropy density at the impact 
parameter b is obtained from the Glauber model 



s{x,y,b) = So 
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where pwN and ps are the densities of wounded nu- 
cleons and binary colUsions respectively. The opti- 
cal Glauber model densities are obtained with Wood- 
Saxon densities for the colliding nuclei pwsi''^) = 
Po/(exp((r-i?,)/a) + l) (po = 0.169fni-3, Ra - 
6.38fni, a = 0.535fni, h — 2.1fni for Au and po — 
0.169fm-3, Ra = 6.48fni, a = 0.535fm, b = 2.2fm for 
Pb) and the nucleon-nucleon cross sections are 42mb and 
62mb at = 200GeV and 2.76TeV respectively. The 
contribution of binary collisions is a = 0.145 at RHIC 
energies [16]. The centrality dependence of the charged 
particle multiplicity in Pb-Pb collisions at the LHC [17] 
is used to fix a = 0.15 at the higher energy. Such a value 
of the parameter a leads to an apparent overprediction of 
the multiplicity in central Pb-Pb collisions [18, 19], but 
this is due to the incorrect value of the reference multi- 
plicity in proton-proton collisions used in Ref. [18]. We 
use a realistic equation of state interpolating between the 
lattice QCD data at high temperatures and an equation 
of state of a gas of hadrons at lower temperatures [20] . 
The use of an equation of state without a soft point is 
essential in reproducing the femtoscopy data [9, 10]. 

The initial time for the hydrodynamic evolution is 
To = 0.6fm/c. The corrections to the pressure from shear 
viscosity 2r]/?>To to the transverse pressure at the initial 
time are of the order of 20%, hence the application of 
viscous hydrodynamics is justified. On the other hand, 
the starting time of the hydrodynamic evolution should 
be defined by the early thermalization mechanisms that 
are not well understood. Moreover, some early flow can 
be generated before the start of the hydrodynamic evo- 
lution. The pre-equilibrium transverse flow in hydro- 
dynamic models has been discussed by several authors 
[13, 21, 22]. The mechanism generating the initial flow is 
not known, and the amount of the pre-equilibrium flow is 
a parameter of the calculation. Ref. [22] gives a general 
result for the transverse flow generated in the initial stage 
of the reaction, which could serve as an upper bound for 
the pre-equilibrium flow generated in the expansion at 
the early stage with boost invariance, starting at zero 
time. 

To estimate the effects of the collective evolution of the 
system before tq and of the choice of the value of tq , we 
compare two different calculation. The first scenario as- 
sumes that the evolution starts at tq = 0.6fm/c with zero 
initial flow and in the second scenario the initial flow is 
taken in the form of the universal flow for the early trans- 
verse acceleration [22]. For a traceless energy- momentum 
tensor the components of the energy momentum tensor 
take a universal form 

at the time tq. 

Starting from tq the evolution is governed by the rel- 
ativistic viscous hydrodynamics with a realistic equation 
of state. The initial energy density and pressure are given 
by the Glauber proflle of the entropy density with cor- 
rections from shear viscosity given by the Navier-Stokes 
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FIG. 1. (Color online) Transverse momentum spectra of tt^ 
in central collisions. The upper lines are for Pb-Pb collisions 
at ^/s = 2.76TeV and the lower lines for Au-Au collisions at 
200GeV. The dashed lines are for standard initial conditions 
and the solid lines are for the evolution starting with the pre- 
equilibrium flow, all calculations in the relativistic viscous hy- 
drodynamics; data for Au-Au collisions are from the PHENIX 
Collaboration [23]. 
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FIG. 2. (Color online)Same as Fig. 1 but for K+ . 

value used as initial conditions for the second order Israel- 
Steward equations. The velocity profile at tq is matched 
to reproduce the ratio T^'/T™ predicted from the uni- 
versality argument (6). The effects of switching to a re- 
alistic equation of state and of the shear viscosity correc- 
tions after tq on the velocity field are small in the core 
of the fireball as compared to using a perfect fluid with 
the ultrarelativistic gas equation of state. 

III. RESULTS 

We set the ratio of the viscosity coefficients to the 
entropy to rj/s = l/Air for the shear viscosity and 
C/s = 0.04 for the bulk viscosity, r^r = rn = The 
bulk viscosity is present only in the hadronic phase [11]. 
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FIG. 3. (Color online) STAR Collaboration data [25] 
(squares) for the HBT radii in central Au-Au collisions at 
= 200GeV compared to the results from viscous hydrody- 
namics with standard (dashed lines) and pre-equilibrium flow 
initial conditions (solid lines). 



FIG. 4. (Color online) ALICE Collaboration data [5] 
(squares) for the HBT radii in central Pb-Pb collisions at 
^/s — 2.76TeV compared to the results from viscous hydro- 
dynamics with standard (dashed lines) and pre-equilibrium 
flow initial conditions (solid lines). 



The statistical emission and resonance decays are per- 
formed using the THERMINATOR program [24], at the 
freeze-out temperatures of 140MeV and 150MeV for the 
calculations without and with with prc-cquilibrium flow 
respectively. The freeze-out with finite shear and bulk 
viscosities in the hadronic phase captures the main char- 
acteristics of the dissipative hadronic rescattering phase 
at the end of the expansion [11]. The entropy density 
at the center of the fireball sq corresponds to a temper- 
ature of 400 and 520MeV at RHIC and LHC energies 
respectively. 

In Figs. 1 and 2 are shown the transverse momen- 
tum spectra of pions and kaons. To compensate for the 
additional transverse flow, the freeze-out temperature is 
increased to 150MeV for calculations including the ini- 
tial flow. The introduction of the initial flow leads to 
some hardening of the spectra as compared to the stan- 
dard initial conditions. However, both calculations are 
close to the experimental data. The transverse momen- 



tum spectra at LHC energies are harder and the relative 
importance of the pre-equilibrium flow is smaller than at 
RHIC. The systems lives longer and most of the trans- 
verse flow is generated in the hydrodynamic expansion of 
an almost perfect fluid. 

For events generated in THERMINATOR the correla- 
tion function is constructed from same-event and mixed- 
event pion pairs [26] . The three dimensional momentum 
correlation functions are fitted using a Gaussian param- 
eterization and the three HBT radii are extracted. The 
dependences of the three radii and of the ratio Rout / Rside 
on the momentum of the pion pair are shown in Fig. 
3 together with the experimental results of the STAR 
Collaboration [25]. The HBT radii are reproduced in 
the viscous hydrodynamic calculation with standard ini- 
tial conditions to within 8-15%, depending on the pair 
momentum. The additional pre-equilibrium flow (solid 
lines) brings the results closer to the data, reducing the 
discrepancy. The satisfactory description of the ratio 
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Rout/ Rside shows that the collective flow generated in 
the model is realistic. The values of the three radii are 
close to the data, which indicates that the size and the 
life-time of the source are well described in the model. 

The source created in Pb-Pb collisions at the LHC is 
larger and shows a stronger collective flow (Fig. 4). The 
viscous hydrodynamic calculation reproduces the mea- 
sured values of the radii to within 14%. The inclusion 
of the initial flow at LHC energies docs not change the 
agreement with the data. The earlier freeze-out for the 
calculation with initial flow leads to a reduction of the 
values of the radii, but the ratio Rout I Rside is not mod- 
ified significantly. We notice that the expansion at the 
LHC is dominated by the hydrodynamic stage, with a 
relatively smaller contribution from the pre-equilibrium 
phase. The strong flow observed in the femtoscopy data 
is quantitatively well described by the model calculation. 
The size of the fireball at the freeze-out increases with 
the increasing multiplicity, as seen in the data. In the 
simulation the freeze-out happens at r ^ 12-13fm/c. 

IV. CONCLUSIONS 

We make hydrodynamic calculations for central heavy- 
ion collisions at = 200GeV (RHIC) and 2.76TcV 
(LHC). The fireball expansion is modelled using relativis- 



tic viscous hydrodynamics with the shear viscosity coef- 
ficient f]/ s = l/47r, and, in the hadronic phase, the bulk 
viscosity of C/s = 0.04 is added. Calculations starting 
without initial flow at tq = 0.6fm/c describe the trans- 
verse pion and kaon spectra at RHIC energies. The HBT 
radii at both colliders are reproduced to within 15%. The 
presence of an additional transverse flow in the initial 
state of the hydrodynamic expansion improves the de- 
scription of the HBT data at RHIC energies. The devi- 
ations of the simulations from the data are of the order 
of the systematic errors quoted by the experimental Col- 
laborations. 

The paper presents the first quantitative description 
of the recently released ALICE Collaboration data on 
the HBT interferometry in Pb-Pb collisions at the high- 
est available energy [5]. The relativistic viscous hydro- 
dynamics gives a satisfactory agreement with the mea- 
surements. The effect of the pre-equilibrium flow on the 
Rout/ Rside ratio is small. 
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